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Quark distribution and spectator functions are estimated in a diquark spectator 
model. The representation of the functions in terms of non-local operators to- 
gether with the rather simple model allow estimates for the yet experimentally 
undetermined functions, like the polarized fragmentation functions or higher twist 
functions. 



More details about the method and the results reported here can be found 
in the long write-up □ by the same authors. 

1 What are the objects of interest ? 

The basic tool in the description of hard processes is the (assumption of) 
factorization into hard and soft physics. The hard parts, i.e., partonic cross- 
sections are calculated using the rules of perturbative QCD, whereas the soft 
parts have either to be taken from experiment or calculated, for instance with 
the help of models. The idea is that the non-perturbative objects containing 
information on the hadronic confinement are universal, i.e., they occur in all 
hard processes in the same form. 

In a field-theoretic description of hard processes the soft, non-perturbative 
physics is contained in (connected) hadronic matrix elements of non-local op- 
erators built from quark and gluon fields. The simplest, but most important 
ones are quark-quark correlation functions. 

1.1 quark-quark correlation functions 

The quark-quark correlation function containing information about how par- 
tons are distributed inside a hadron with momentum P and spira. S in the 
initial state of the process can be written for each quark flavor as Ha 

$«(p,P,5) =J2j-0f (P,-5|^.(0)|X)(X|^(x)|P,5). (1) 
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The information on how a quark, with a specific flavor, fragments into a hadron 
with momentum Ph and spin Sh (plus anything else; denoted by X) is encoded 
in the correlation function □ 



x •* 



(2tt)< 



e lk - x (0\^(x)\P h ,S h ;X)(P h ,S h ;X\ijJQ)\0). 



(2) 



To be specific we consider the case were the hadrons are protons. 
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1.2 distribution junctions 

In a particular hard process only certain Dirac projections of the quark-quark 
correlation functions appear and the non-locality is restricted by the integra- 
tion over quark momenta to light-like directions. 
The Dirac projections of <&ij(p, P, S) 



$[ r l( x ) = I J dp- d 2 p T Tr (*r) 



(3) 

cP+ 



define the usual distribution functions. The functions /i, g\ and h\ obtained 
by 

t> h+] (x) = fi(x) (4) 

$[7 + 75]( a .) _ \ gi ( x ) (5) 

^ a+ ^][x) - S%hi(x) (6) 

are leading in an 1/Q expansion, where Q is the hard scale typical to the pro- 
cess. A is the helicity of the hadron and St the transverse part of the spin 
vector. Leading functions have an intuitive probabilistic interpretation: 



• /i (x) gives the probability of finding a quark with light-cone momentum 
fraction x in the + direction (and any transverse momentum). 

• gi (x) is a chirality distribution: in a hadron that is in a positive helicity 
eigenstate, it measures the probability of finding a right-handed quark 
with light-cone momentum fraction x minus the the probability of finding 
a left-handed quark with the same light-cone momentum fraction. 

• h\ (x) is a transverse spin distribution: in a transversely polarized hadron, 
it measures the probability of finding quarks with light-cone momentum 
fraction x polarized along the direction of the polarization of the hadron 
minus the probability of finding quarks oppositely polarized. 

The sublcading ('higher-twist') functions have no intuitive partonic interpre- 
tation. Nevertheless, they are well defined as projections of the quark-quark 
correlation functions. The pre-factor M/P + behaving like l/Q in a hard pro- 
cess signals the sub-leading nature of the corresponding distribution functions 

= ^e(x) (7) 
$^(x) = ^SZg T (x) (8) 
<S> [ta+ ~~< 5] (x) = M.\h L {x). (9) 

1.3 fragmentation functions 

Analogously, the fragmentation functions are defined by the Dirac projections 



A^z) = ±-J dk + Gpk-rTr (AT) 



(10) 



k-= P -/z 



where the leading ones are obtained by 

A^(z) = D^z) (11) 
A^^(z) = A fc Gi(*) (12) 
A [^-75] (z) = SfcH^z) (13) 

and sublcading 'twist' fragmentation functions by 

AW(z) = ^E(z) (14) 
Mi 
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A^(z) = L±S%G T (z) (15) 

A^- + ^(z) = MH\ h H L (z). (16) 

Note the pre-factor Mh/P^ (scaling like l/Q) accompanying the latter 
functions. 

There are 3 additional subleading fragmentation functions: Dt{z), E^{z) 
and H(z), which are so-called "T-odd" functions Since we will describe the 
states of the produced hadrons by free spinors, our model leads to identical 
vanishing T-odd functions. 



2 Spectator model 

2.1 Strategy to obtain estimates 

The purpose of our investigation is to obtain estimates for those distribution 
and fragmentation functions which are experimentally poorly (or not at all) 
known at present. These are, in particular, the transverse spin distribution 
hi, the polarized fragmentation functions Gi, Hi and the subleading twist 
functions e(x), gr(x), hhix) and E(z), Gt(z), Hl(z). 

To this end we employ a rather simple spectator model with only a few 
parameters. After fixing the parameters by phenomenological constraints we 
check that the gross features of the experimentally well-known distribution 
functions fiix) and gi{x) and the fragmentation function Di(z) are satisfac- 
torily reproduced. 

The advantages of the model are its Lorentz covariance and its easy exten- 
sibility to polarized and subleading twist functions. We end up with analytic 
expressions for all the desired functions. 



2.2 ingredients of the model 

The ingredients of the model are exemplified for the distribution functions [| 

• The basic idea of the spectator model is to assign a definite mass to the 
intermediate states occurring in the quark-quark correlation functions 
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(P-p) 2 =M R 2 



(k-P h ) 2 =M R 2 



• The quantum numbers of the intermediate state are determined by the 
action of the quark field operator on the hadronic state \P, S), i.e. they 
are the quantum numbers of a diquark system 
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The matrix element appearing in the RHS of (Rl) is given by 



{x,mo)\p,s) 



T s kl U t (P,S) 



ik 



(17) 



in the case of a scalar diquark, or by 



(x*\mo)\p,s) = 



y 1 — m 



T°f Ut(P,S) 



(18) 



for a vector diquark. The matrix elements consist of a nucleon-quark- 
diquark vertex Y(N), the Dirac spinor for the nucleon Ui(P, S), a quark 
propagator for the untruncated quark line and a polarization vector e 
in the case of an axial vector diquark. 



• For the nucleon- quark- diquark vertex we assume the following Dirac struc- 
tures []: 



quark 




diquark 

(singlet) 



quark 



diquark 

(1 

(triplet) 



T s (iV) = 1 g s (p 2 



T afl (N) 



la( P 2 ) 



v/3 



pp. 

a7 



75 [I 1 * + T7 ( 20 ) 



(19) 



c a special case of the most general formL 



• The functions g s / a {p 2 ) are form factors that take into accaunt the com- 
posite structure of the nucleon and the diquark spectator 0. We use the 
same form factors for scalar and axial vector diquark: 

9sfa( P 2 )-N^-^. (21) 

• The flavor coupling of the proton wave function from a scalar diquark 
(Sq) and an (axial)vector diquark with isospin or 1 (Aq or A\, respec- 
tively) 

\p) = -L\uS ) + -^=\uA )--L\dA 1 ), (22) 
leads to the flavor relations 

fx = \fi + \ft ( 23 ) 
ft = ft (24) 

and similarly for the other functions. The coupling of the spin has already 
been included in the vertices. 

Putting all ingredients together analytic expression for the quark-quark corre- 
lation functions are obtained 

<f>n (n P q\ N " ^~2P-p + NP-Ml) 

\P,r,D) 2(2?r) 3 | p 2_ A 2|2a 

x(J + m){^+M)(l + a Rl5 $)(^ + m) (25) 

from which the distribution functions can be easily projected. The calculation 
of the quark-quark correlation function for the fragmentation, A^-(fc, P^, Sh), 
proceeds along the same line. 



2.3 fixing the parameters 

The parameters of the model are fixed as follows: the power in the denominator 
of the form factor, a = 2, is chosen to reproduce the Drell-Yan-West relation for 
large x. The mass difference M a — M s is motivated by the TV— A mass difference 
(with group theoretical factors properly accounted for) and the values for A 
and M s reproduce the experimental value for the axial charge, qa- 
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parameter 


fixed by: 


a = 2 


Drell-Yan-West relation 


M a -M„ = 200 MeV 


N — A mass difference 


M s = 600 MeV; A = 0.5 GeV 


g A = 1.28 


N 


number sum rules 




fdxf? = 2;fdxf? = l 



3 Numerical results 

Having fixed the parameters of the model we can present the numerical results 
for the distribution and fragmentation functions* Analytical expressions for 
the functions can be found in the long write-up 111. Note that the model does 
not provide a scale dependence; but it is expected to describe physics at a low 
"hadronic" scale. Thus we compare to distributions found in the literature at 
scales as low as available. 

• The experimentally best known function is the distribution fx(x). We 
compare our results with the parametrization from Gliick, Reya and 
Vogt El at the low scale /i| = 0.23 GeV 2 . Note that the first mo- 
ments, J fi(x)dx are exactly the same (normalization condition), a fact 
not immediately apparent from the diagram, since we plot the combina- 
tion x * f\ (x) . Our distribution is narrower due to the non- inclusion of 
gluons and anti-quarks (which - if included in our model - would have a 
broadening effect). Thus we refrain from fine-tuning parameters to ob- 
tain a better agreement with GRV — and are satisfied with agreement 
in the gross features. 
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Less well experimentally determined is the helicity distribution gi(x). 
Our predictions are again smaller than the distribution taken from the 
literature E21 Signs and positions of extrema are in reasonable agreement. 




x g,(x) 



it - quark 
d - quark 



x Ah/jmW )[GRSV] 
xAdJx.jlJUGRSV] 




Yet experimentally undetermined is the transverse spin distribution hi (x) . 

Within our model we obtain an 
xh,(x) » estimate for this function which 




is numerically close to the re- 
sult for the helicity distribution 
gi(x), but not identical. 



As well unknown are the subleading 'twist' functions like e(x) and the 
combination 172(2;) = gr(x) — gi{x), for which we display our estimates 
below. 




On the side of the fragmentation functions only the spin-independent 
information contained in D\(z) has been experimentally investigated. 



Unfortunately, no parametrizations at a low "hadronic" scale are avail- 
able in the literature. Thus we compare our findings with data from the 
EMC collaboration^ at 20 GeV 2 ; keeping the mentioned shortcoming in 
mind. We plot for the data the combination Z?"^ p — D^~ yp . The latter 
term being predicted zero in our model (no anti-quarks included !) this 
"valence" combination seems to be best suited for a fair comparison. 

A second difficulty is caused by the fact that for fragmentation functions 
we do not have a number sum-rule to fix the normalization. Thus we 
normalize our findings such that they reproduce the second moment, 
J z D\(z)dz (a scale invariant quantity) as obtained from the EMC data. 

We avoid the normalization problem for the (longitudinal and transverse) 
spin-dependent fragmentation functions, G\(z) and Hi(z), by showing 
their ratios with D\{z). 
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• The subleading 'twist' fragmentation functions have an even worse chance 
to be experimentally investigated in the near future than the leading spin- 
dependent ones. We don't display numerical results for them, although 
they are easily obtainable in the context of our model. All necessary 
information to write down their analytic expressions is contained in our 
article EL 

4 Conclusions 

We presented estimates for distribution and fragmentation functions obtained 
in the context of a spectator model using representations of the functions in 
terms of non-local operators. 

By comparing our expressions to available parametrizations at low "ha- 
dronic" scales and to some experimental data we find that we can obtain 
reasonable qualitative agreement for fi(x), gi(x) and D\(z). These findings 
give confidence that the estimates obtained for the "terra incognita" functions 
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(transverse spin distributions, longitudinal and transverse spin fragmentations 
and subleading functions) provide a reasonable estimate of the order of mag- 
nitude of the functions and their (large) x behavior despite the simplicity of 
the model. 

Moreover, the model predictions fulfill not only trivial positivity con- 
straints like < fi(x) and |/ii(x)| < fi(x), but also for instance the 
Soffer inequalityE3, 2|/ii(x)| < (fi(x) + gi(x)); the latter in our case becoming 
an equality []. 
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